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Abstract 

Recently, models proposing superheavy particles X as source of ultrahigh energy 
cosmic rays have attracted some interest. The X-particles are either metastable 
qq " relic particles from the early Universe or are released by topological defects. In 

these models, the detected air-showers are produced by primaries originating from 
the fragmentation of the X-particles. We present the fragmentation spectrum of 
superheavy particles calculated in SUSY-QCD. Then we discuss the status of the 
lightest supersymmetric particle as possible ultrahigh energy primary. 
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1 Introduction 



Cosmic raysB (CR) are observed in a wide energy range, starting from subGeV energies 
up to 3 ■ 10 20 eV. Apart from the highest energies, these particles are accelerated in 
our Galaxy, most probably by shocks produced by SNII explosions. There is no universal 
definition of Ultra High Energy Cosmic Rays (UHECR). We will use this term for energies 
E > 10 19 eV, where a new, natter component appears in the CR spectrum (Fig. 1). The 
highest energies detected so far are 2 — 3 • 10 20 eV. 
Ph| It is natural to think that the UHE component has an extragalactic origin, since the 

galactic magnetic field cannot confine particles of these energies. Moreover, the accel- 
eration of protons or nuclei up to 2 — 3 • 10 20 eV is difficult to explain with the known 
astrophysical galactic sources. The most prominent signature of extragalactic UHECR 
is the so called Greisen-Zatsepin-Kuzmin (GZK) cutoff B: The energy losses of protons, 
nuclei and photons sharply increase at -Egzk ~ 3 • 10 19 eV, reducing the mean free path 
length I of these primaries to less than 50 Mpc or so. Thereby, the spectrum should be- 
come stepper above -Egzk for any source with a distribution homogeneous on scales larger 
than I. There is another argument which disfavours the standard astrophysical sources: 
At energies E ~ 10 20 eV, the arrival direction of the primaries (which is known within 
several degrees) should point towards their site of origin. But no source of UHECR like, 
e.g., active galactic nuclei has been found within 50 Mpc in the direction of these events. 

An elegant solution to the above problems are top-down models: In contrast to the 
standard sources, the primaries are not accelerated but are the fragmentation products 
of some decaying superheavy particle X. For X-particles with mass nix 3> -^gzk, the 
acceleration problem is solved trivially. Moreover, these sources also evade detection by 
normal astronomical methods. 

This contribution is organized as follows: In Section 2, the in our view two most 
promising top-down models are presented. In Section 3, the spectrum of hadrons pro- 
duced in the decay of supermassive particles is calculated in supersymmetric QCD. As an 
application, the CR fluxes produced by cosmic necklaces are shown in Section 4. Finally, 
the status of the LSP as UHE primary is discussed in Section 5. 
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Figure 1: Compilation of the CR spectrum from AGASA, Fly's Eye, Haverah Park and Yakutsk. A two 

component fit to the spectrum is also shown. 

2 Top-down models as new sources for UHECR 

Three different possibilities for top-down models have been discussed in the literature: 

1. Primordial black holes: During the final stage of their evaporation, they emit high 
energy particles. However, when the resulting spectra of cosmic rays are combined 
with various observational bounds on the mass fraction of the universe in black 
holes, one finds that the UHE CR flux from black holes is well below the observed 
flux. 

2. Topological defects I such as superconducting strings, monopoles, and monopoles 
connected by strings: Here we will concentrate on cosmic necklaces, since this model 
seems to provide the largest UHE particle flux for fixed density of electromagnetic 
cascade radiation. 

Cosmic necklaces are hybrid defects consisting of monopoles connected by a string. 
These defects are produced by the symmetry breaking G — > H xU(l) — > H x Z2. In 
the first phase transition at scale rj m , monopoles are produced. At the second phase 
transition, at scale r] s < r] m , each monopole gets attached to two strings. The basic 
parameter for the evolution of necklaces is the ratio r = m/(/i<i) of the monopole 
mass m and the mass of the string between two monopoles, /i<i, where fi ~ rj 2 is 
the mass density of the string and d the distance between two monopoles. Strings 
loose their energy and can contract due to gravitational radiation. As a result, all 
monopoles annihilate in the end producing superheavy Higgs, gauge bosons and 
their supersymmetric partners which we call collectively A-particles. The rate of 
A-particle production can be estimated as 

dn x ^ r 2 \x 
dt ~ t 3 m x ' 

The flux of UHECR is determined mainly by two parameters, r 2 fi and mx, which 
values must be of order 10 27 GeV 2 and 10 14 GeV, respectively, to have the flux close 
to the observed one. For a more complete discussion see Ref. [4]. 



3. Superheavy, metastable relic particles X: They constitute (part of) the cold dark 
matter (CDM) and, consequently, their abundance in the galactic halo is enhanced 
by a factor ~ 5 x 10 4 above their extragalactic abundance. Therefore, the proton 
and photon flux is dominated by the halo component and the GZK-cutoff is avoided 
as was first pointed out in Ref. [5]. 

The necessary lifetimes, r > 10 17 s, and mass ranges, m>x > 10 13 GeV, of the 
X-particle arise quite naturally in several extensions of the standard model. The 
X-particle could be protected by some sort of i?-parity which is extremely weakly 
broken by wormholei or instantoni effects, or could belong to the hidden or messen- 
ger sector of SUSY mo dels!. The most promising production mechanism proposed 
so far is the enhancement of vacuum fluctuations of the X field during inflation 
by gravitational interactions!. For 10 12 GeV < nix < 10 13 GeV, this mechanism 
results in a V-particle abundance close to the critical one. 

A common signature! of all top-down models is the high photon/proton ratio 7/p. 
Since in the fragmentation process much more mesons than baryons are produced, the 
ratio 7/p should be > 1 at the highest energies. (The exact value depends on the frag- 
mentation model and on the poorly known absorption length of UHE photons.) If the 
UHE primaries originate from the decay of some superheavy CDM particles, an additional 
signature is their anisotropyl'0 reflecting the non-homogeneous distribution of the CDM. 

3 Fragmentation spectrum of hadrons in SUSY-QCD 

The spectra of hadrons produced in deep-inelastic scattering and e + e~ annihilation are 
formed due to QCD cascading of the partons. In the Leading Logarithmic Approxi- 
mation (LLA) which takes into account \n(Q 2 ) terms this cascade is described by the 
Gribov-Lipatov-Altarelli-Parisi-Dokshitzer (GLAPD) equation0. This approximation is 
not valid for Bjorken x< 1, when colour coherence effects become important. A bet- 
ter approximation is the Modified LLA (MLLA), which takes into account both ln(Q 2 ) 
and ln(x) terms as well as angular ordering. In this section, we consider the limiting 
spectrum @ that has been obtained as an approximate solution to the MLLA evolution 
equations for small x. In fact, it describes well also the experimental data at x ~ 1. 

At large energies > 1 TeV the production of supersymmetric particles might 
substantially change the QCD spectra. Apart from future experiments at LHC, super- 
symmetry (SUSY) might strongly reveal itself in the decays of superheavy particles: As 
long as the virtuality Q 2 remains much larger than the SUSY scale M| USY , in the particle 
cascade initiated by the decay of a V-particle not only usual particles but also their super- 
symmetric partners participate. As we will see, the fragmentation spectrum of hadrons 
changes considerably going from QCD to SUSY-QCD lid. To obtain reliable results for 
UHECR fluxes, it is therefore necessary to calculate the spectra within SUSY-QCD. 

We consider first the SUSY-QCD cascade in LLA. Neglecting terms proportional to 
a{Q 2 ) and keeping terms with a s (Q 2 ) ln(Q 2 ), the GLAPD equation can be written aslll 

-57^,0 = £/ -<S> C a (z)DE(x/z,0-Y< dzz<l> c A (zM(x,0, (2) 
ot, q Jo z c 

where &a( z ) * s the splitting functional characterizing the decay A — > B + C. Here 
D\ is the distribution of partons B inside the parton A dressed by QCD interactions 



with coupling constant a(k\), where k± is the transverse momentum and x = k\\/k™ x is 
the longitudinal momentum fraction of the parton B. The variable £ characterizes the 
maximum value of k 2 L available in the considered process (k\ < Q 2 ), 

t(cfl\ - f Q2 dk ^ a ^ (*\ 

with A ~ 0.25 GeV as phenomenological parameter. 

The supersymmetrization of Eq. (^) is simple: each parton A should be substituted by 
the supermultiplet which contains A and its superpartner A. We are mainly interested in 
the spectrum at small x, where gluons strongly dominate the other partons. Therefore, it 
is a reasonable approximation to consider only two partons, namely gluons g and gluinos 
A, in the tree diagrams. 

After performing a Mellin transformation, 

DiU,Z)= f'dxx^D^x^), (4) 
Jo 

Eq. (|2|) can be rewritten in matrix form, 

-^D(j,£) = H(j)D(j,0, (5) 
where we have chosen as basis (g, A) and 




u 9 (j) = j Q dz [[z^ - z) & g (z) - ^ g (z)\ (7) 

u x {j) = £ dz (z^ 1 - l) $ x (z) . (8) 
After diagonalization of H(j), the eigenvalues of H are 
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v± = -[v g + u x ± (u g -u x Y + A^{ ). (9) 
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In the limit u = j — 1 —>■ 0, the leading term v + is given by 



v+= — - - a + 0{u) (10) 
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with a = yiV c = 11. (N c = 3 is the number of colours.) 

Dokshitzer and TroyanEl were able to express the limiting spectrum D\i m which is a 
MLLA result as a function of the parameter a calculated in LLA and b, the constant of 
evolution of a s (k±) in one-loop approximation, k^dasik 2 ^) / dk\ = — 6af/(47r). Properly 
normalized, the function Z} lim (Z,F) = xD\ im (x,Y) gives a~ l da/dl in the case of e + e~ 
annihilation and the decay spectrum of X particles, where / = ln(l/x), Y = ln[-y/s/(2A)] 
and \/i is the cm. energy of an e + e~ pair or the mass nix of the superheavy decaying 
particle. It is given by 

^ /, *rs ^ 4Cp„,„, r /2 dr Rn ( b sinha Y\ B r , . , „, 

DUl,Y)=K lim —^T(B) e- Ba [ — I B (y) 11 

o j-tt/2 7r \oiv c a y I 



with a = 11 and b = &susy = 9 — n/ = 3in SUSY-QCD. Furthermore, a = a + ir, 
a = arctanh(2C - 1), C = 1 - l/Y and C F = (N? - 1)/(2N C ) = 4/3. Finally, I B is the 
modified Bessel function of order B = a/b, and argument 

y{r) =[ \ si ^ ha [cosh a + (1 - 2() sinha] Y) . (12) 

In Fig. 2, the SUSY-QCD limiting spectrum is shown in comparison with the QCD limiting 
spectrum. The maxima of the SUSY spectra are shifted to the right, and, since they are 
also narrower than the QCD spectra, the SUSY maxima are dramatically higher. 
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Fig. 2 (left): Limiting spectrum D Xim (l,Y) for SUSY-QCD (solid lines) and QCD (dashed lines). The 
QCD spectrum is scaled up by a factor 30. Both cases for mx = 10 12 GeV (bottom), mx = 10 13 GeV 
(middle) and m x = 10 14 GeV (top). 

Fig. 3 (right): Predicted fluxes from cosmic necklaces with r 2 /i = 5 • 10 27 GeV for mx = 10 14 GeV. 



4 Application: UHECR fluxes from necklaces 

We present now the application of the SUSY-QCD limiting spectrum to the calculation 
of UHECR spectra generated by the decay of superheavy particles^. 

Let us discuss first the problem of the normalization of spectrum. The normalization 
constant Ku m cannot be calculated theoretically and is normally found from comparison 
with experiment. Since the spectrum changes dramatically going from QCD to SUSY- 
QCD, we cannot use this value of K\ im . Instead we use as normalization condition 

/ dx xD i]5m (x, Y) = 2fi , (13) 
J o 

where % runs through iV (all nucleons) and ir^ and 7r° (charged and neutral pions), while 
fi is the fraction of energy carried by the hadron %. Note that the main contribution 
to the integral in Eq. ( |I~3"D comes from large values x ~ 1, where the limiting spectrum 
might have large uncertainties. However, this is, in our opinion, the most physical way of 
normalization. The numerical values of f\ are unknown at large s. One can assume that 
f n ~ 1 — /lsP) where /lsp is the energy fraction taken away by the lightest supersymmetric 
particle (LSP). According to a simplified Monte-Carlo simulation Ej, /lsp ~ 0.4. For the 
ratio /jv//tt we use ~ 0.05 inspired by Z° decay. 



Let us assume that the decay rate of X-particles fix in the extragalactic space does 
not depend on distance and time. Then taking into account the energy losses of UHE 
protons and the absorption of UHE photons due to pair production (7 + 7 — > e + + e - ) on 
the radio and microwave background, the diffuse flux of UHE protons and antiprotons is 



p+p 



1 h x 



2n mx Jo 



dt g D N (x g ,Y) 



dEg(E, tg 

dE 



(14) 



where E g (E,t g ) is the energy at generation time t g of a proton which has at present the 
energy E and x g = 2E g /mx- Denoting the proton energy looses on microwave radiation 
by dE/dt = b(E, z), dE g /dE is given by 



dE g (E,z g 
dE 



1 + z g ) exp 



E=E g (z). 



(15) 



where Ho is the Hubble constant and z the redshift. The diffuse spectrum of UHE photons 
can be calculated as 



I {E) = \ [E) J- f 1 — D 7r o(x,Y) 

7T m X J2E/m x X 



(16) 



where A 7 is the absorption lengthlll of a photon. Finally, the diffuse neutrino flux is gi 
by 



iven 



;i7) 



where 1 + z max (E) « mx/(4i?). The neutrino flux depends generally on the evolution of 
the sources, 



ft \ ■ 

hx(t) = h x {to) ij) 



;is) 



In Fig. 3, the spectra of UHE protons, photons, and neutrinos are shown together 
with experimental data for the model of cosmic necklaces with r 2 fi = 5 ■ 10 27 GeV 2 and 
mx = 1 • 10 14 GeV. The proton flux is suppressed at the highest energies as compared 
with the calculations of Ref. [4], where the Gaussian SUSY-QCD spectrum was used. 



5 The LSP as UHE primary 

We have noted in Section 3 that usual particles as well as their super symmetric partners 
participate in a particle cascade as long as their virtuality Q 2 remains larger than the 
SUSY scale M| USY . However, when Q 2 reaches M| USY , the supersymmetric particles 
stop branching and decay to the LSP. Thereby, UHE LSP will be produced in top-down 
mo dels@ (assuming i?-parity is unbroken). Since the LSP cannot be effectively accelerated 
in standard astrophysical sources, the detection of UHE LSP in extensive air shower (EAS) 
experiments would not only show that SUSY is realized by Nature but would be also a 
clear signature for top-down models. 



5.1 Gluino as LSP 



Hadronic bound states of the gluino g (which we call generically g-hadrons G) were 
discussed already in the 80's as primary of the UHECR@. Recently, this idea was revived 
by Farrar et al. proposing as UHE primary the gluebarino S = guds. A simple model 
that leads to the MSSM with gluinos as LSP was presented by S. Raby0. 

There are several arguments against a very light gluino with 0(m§) = 1 — 10 GeV: 
First, the running of a s gives the bound m- g > 6.3 GeV. Second, two searches for g- 
hadrons at Fermilab had negative results HI. Third, g-hadrons are produced by CR in 
the earth atmosphere. In the case that the lightest gluebarino is the state guud and 
stable, the light gluino as LSP is excluded by the search for heavy hydrogen c3. If the 
lightest gluebarino is neutral, as considered by Farrar et al., this argument HI still work 
if S forms a bound state with the nuclei. Thus, in our view, a very light gluino is 
disfavoured. If the gluino is heavy enough (according Ref. [20] m- g > 50 GeV), it will 
produce the standard missing energy signal in accelerator experiments. Therefore, the 
range 50 GeV < m 5 < 154 GeV is also excluded. 

The interactions of UHE g-hadrons were already considered for the case of a glue- 
ballino. Two values determine the interaction of a UHE g-hadron with a nucleon. The 
first one is its total GN-cross-section which should be proportional to its size. However, 
different estimates have been made in the literature: a ~ a s (mg)/fi 2 ~ 1 mb@, where u 
is the reduced mass of G, a ~ A~^ CD ~ 10 mb0, and a ~ (1 — 10) a n N ~ (3 — 30) mbH. 
The second important value is the average energy transfer (y) per scattering. For the 
production of EAS in the atmosphere only interactions with large y are effective: While a 
very light g-hadron would interact like a nucleon, a heavy one behaves in the atmosphere 
like a penetrating particle and should be distinguishable from a proton 011. 



5.2 Neutralino as LSP 

In most SUSY models, the neutralino is the LSP and, in contrast to the gravitino and 
gluino, it is also a viable DM candidate. Let us consider the interactions of the neutralino \ 
relevant for their detection^. Mainly two processes are important for its interaction with 
matter, namely the neutralino-nucleon scattering \ + N — > all and resonant production 
of selectron off electrons x + e — > & ~ > an - The first process is based on the resonant 
subprocess x + g — > g all and on neutralino-gluon scattering. The latter subprocess is 
important, because for high energies, and consequently for small scaling variable x, the 
gluon content of the nucleon increases fast. 

The cross-sections of all subprocesses start to grow at energies s 3> m|. The rise with s 
is caused by the decrease of x m i n = mhjs and x m i n = (m,- + m g ) 2 /s, and the corresponding 
increase of the number of partons with sufficient momentum in the nucleon. If squarks do 
not decay mainly into neutralino, i.e. T tot 3> T (<?£,# — > <? + x), neutralino-gluon scattering 
gives the dominant contribution to the total cross-section. At energies s m 10 10 (GeV) 2 
or E x ~ 5 • 10 18 eV, the neutralino-nucleon cross-section is about 10 -35 — 10~ 34 cm 2 (for 
irtq ~ 900 GeV), i.e. slightly lower than the neutrino- nucleon cross-section. 

Let us consider now x+e — > e — > all which is similar to the Glashow resonant scattering 
z/ e + e— >IU— > n + Vp. The resonant energy of the neutralino is 



E x = M?/(2m e ) = 9.8 ■ 10 8 (M e -/10 3 GeV) 2 GeV 



(19) 



and the cross-section is given by the usual Breit-Wigner formula. The resonant events 
are produced as a narrow peak at E x and give an unique signature for the neutralino. 



6 Conclusions 

Superheavy, metastable relic particle and cosmic necklaces are the two most promising 
sources for top-down models. Their signature is the high photon/proton ratio and the 
LSP as UHE primary. The UHECR flux produced by relic particles has additionally a 
small galactic anisotropy and, as most prominent signature, no GZK-cutoff. For a reliable 
calculation of UHECR fluxes the knowledge of the fragmentation spectrum of superheavy 
particles is necessary. We presented the limiting spectrum in SUSY-QCD that differs 
drastically from the QCD spectrum. 
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